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Abstract 



The Euler equations associated with difteomorphism groups have received much recent study 
because of their Hnks with fluid dynamics, computer vision, and mechanics. In this paper, we 
consider the dynamics of N point particles or 'blobs' moving under the action of the Euler 
equations associated with the group of diffeomorphisms of the plane in a variety of different 
metrics. The 2 body problem is always integrable, and we analyze its phase portrait under 
different metrics. In particular, we show that 2-body capturing orbits (in which the distances 
between the particles tend to as t — *■ oo) can occur when the kernel is sufficiently smooth 
and the relative initial velocity of the particles is sufhciently large. We compute the dynamics 
of these 'dipoles' with respect to other test particles, and supplement the calculations with 
simulations for larger N that illustrate the different regimes. 

1 Introduction 

The Euler equations for planar diffeomorphisms are 



where m denotes differentiation with respect to time, u{x,t) {u,x e M", i e K.) is a velocity 
field, and m{x, t) its associated momentum. The velocity u and momentum m are related by 



where A is an elliptic operator (e.g. {\ — V^)*^) called the inertia operator. 

There has been much recent interest in the Euler equations because they arise in several 
different fields, including computer vision and fluid dynamics. In computer vision they appear 
in two and three dimensions in the field of image registration, such as in the averaged template 
matching equations ^1 and the geodesic interpolating clamped-plate spline , while in 
fiuid dynamics they are the limiting case of the shallow water wave Camassa-Holm equations 
in one and two dimensions |7j, and the motion of an ideal incompressible fiuid corresponds to 
geodesies on the group of volume-preserving diffeomorphisms ini 

Euler equations such as Eq. have a natural geometric origin, which we sketch here 
although these details are not needed in the paper. Let © be a Lie or diffeomorphism group 
with Lie algebra g. Let © be equipped with a left- or right-invariant metric that restricts to a 
metric ((,)) on q. Typically, this metric is defined by a linear inertia operator A:q^ Q* via 



m -\- u ■ Vm + Vii"^ • m -|- m(div u) = 0, 



(1) 



m — Au 



(2) 



{{u,v)) := {u,Av). 



(3) 
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The geodesic equation on T(& can be reduced to give a non-canonical Hamiltonian evolution 
equation on g called the Euler-Poincare equation, or transferred (via the Legendre transform, 
0^0*, which in this case is u — > m := Au) to a Lie-Poisson system on g* called the Euler 
equations 

TO = ±ad;^-i^TO, (4) 

where the sign is + for left- and — for right-invariant metrics. The variable m e g* is called the 
momentum and u € Q the velocity. 

The most famous Euler equation on a Lie group is the equation of motion of a free rigid 
body (© — 50(3), m — body angular momentum). An infinite-dimensional example is the 
Landau-Lifshitz equation on the loop group (S = C°°(M, 50(3)). Examples on diffeomorphism 
groups include the Euler fluid equations on the group = Diffvoi(K") of volume-preserving 
diffeomorphisms with respect to the metric (to = vorticity) 5 , the Camassa-Holm equation 
((S = Diff(§i), metric) H, and the second-grade fluid equations (25 = Diffvoi(M"), 
metric) 

There has been recent interest in the use of diffeomorphic deformations ('warps') to align 
groups of medical images. This alignment allows the range of variation in the appearance of 
biological structures to be measured, providing an approach to automated disease diagnosis. 
One approach to this problem, known as Computational Anatomy, is to introduce group actions 
as deformable templates that are warped via the actions of a group onto other images [HI ■ This 
work has been fundamental to a large amount of research on aligning images through landmark 
matching, where corresponding points arc defined on a set of images, and diffeomorphic warps 
used to align them. The corresponding group for image analysis is the full diffeomorphism 
group, not the volume-preserving subgroup. Under a right-invariant Riemannian metric, it 
can be shown that the geodesies of the motion of a set of landmarks can be computed as an 
optimisation problem. See |12l I13j for an overview. 

For 6 ~ Diff(R") the Euler equations are given by Eqs. (|1I2|I (see [HllEl for further details). 
The inverse of the inertia operator A is given by convolution with the Green's function G of A, 
i.e., u — G *m, where * denotes convolution and AG{x, x') = 5{x — x') for x, x' € R". We shall 
only consider rotationally invariant and diagonal A\ in this case G{x,x') — G{\\x — x'\\) for a 
scalar function G, which we call the kernel of G. Examples in the literature include Gaussian 
in] and Besscl function kernels ^] associated with the metric. In this paper, we consider a 
family of kernels associated with metrics for various fc, which includes the Gaussian kernel 
in the limit as fc ^ oo. 

A striking feature of Euler equations on diffeomorphism groups is that they admit (formally, 
at least) exact solutions in which the momentum is concentrated at a finite set of points that 
we call particles. For fluid equations these are point vortices, which are widely studied both in 
their own right and as a means of approximating the evolution of smooth or other vorticities. 
At flrst sight it is remarkable that a PDE should have such a flnite-dimensional reduction, 
unrelated to integrability. (Such a reduction, however, is not unique to Euler equations. On the 
other hand, there are many particle-based numerical methods, collectively known as Smoothed 
Particle Hydrodynamics ^3], which are not exact solutions of the PDE they approximate.) 

For the 2D and 3D Euler fluid equations, convergence of the point vortex solutions to solu- 
tions for smooth initial data has been established ^ . The speed of convergence can be improved 
by smoothing out the point vortices to vortex blobs even though the (e.g., Gaussian) blobs 
are no longer an exact solution of the Euler equations. Instead, their evolution can be regarded 
as that of delta- functions under a slightly different inertia operator. Altering the shape of the 
blob (which corresponds to the kernel) is equivalent to altering the metric. In this way we are 
lead to consider a wide class of metrics. For an introduction to many-body problems, particu- 
larly arising in Newtonian and relativistic mechanics, see Point vortices defined over Dirac 
delta- functions, and the regularised version, vortex blob methods, are introduced in j^]; see [TU| 
for a review of their use. 

The methods used for image registration and template matching of images in |13l IT^ , and 
references therein, correspond to computing the motion of particles in Diff(R^). Various metrics 
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have been used, including Gaussians, H^, and H'^. DifFcomorphic image registration can be 
considered as a boundary value problem, since the start and end points of the particles are 
defined by the images under consideration. However, it can equivalcntly be treated as an initial 
value problem under the definition of the starting point and initial velocity of the particles |17l 
I12j . This paper provides an initial study of the dynamics of particles for the group Diff(IR^) of 
all diffeomorphisms of the plane, as are used in the image registration framework. Following the 
approach taken in fluid dynamics, we consider a set of N particles moving under the action of 
a kernel G. 

We begin the paper by explicitly computing the Hamiltonian for the case of N particles, and 
computing the equations of motion in terms of the Green's function. Following this, in section|31 
we use delta-function particles and Helmholtz-style inertia operators to find the kernels Gk 
corresponding to H'^ metrics. The limit of these, H°°, corresponds to a Gaussian kernel for 
delta-function particles, where the inertia operator tends to exp(— e^V^). We show how even 
when blobs are required, rather than particles, we can still consider delta-function particles 
through a change in the inertia operator. 

We then shift our attention to the case N = 2, and examine the simplest possible interaction 
of particles. This two body problem is always integrable. We derive the reduced Hamiltonian 
and consider how the particle dynamics change for different metrics. The most striking feature of 
the dynamics is the phenomenon of particle capture, in which particles cohere or stick together. 
In a typical capture orbit, the distance between two particles is ©(e""*) as i ^ oo, but there is 
no actual collision. We identify the regimes of capture and (more traditional) scattering orbits. 
For different metrics, we also show the phase portraits of the dynamics, and the appearance of 
typical scattering orbits. The analysis shows that the phase portrait depends sensitively on the 
behaviour of G close to r = 0. Only when the kernel G is smooth at r = and the relative initial 
velocities of the particles are sufficiently large will the particles capture each other, otherwise 
they scatter. For the cases k < 2 the Green's function is not sufficiently smooth to typically 
allow capture orbits. 

Following this, in section|Slwe consider the case fc > 2, so that the kernel is sufficiently smooth 
to allow capture orbits for 2 particles, and show how the dynamics of a set of N particles varies 
as the particles approach each other. We consider the behaviour of 'dipoles', as we term pairs of 
particles that reach the limiting state of capture, and show examples of typical orbits for a group 
of A'' = 20 particles over a number of timesteps. We conclude with a summary of the differences 
between the particle dynamics in Diff(]R^) with an H'' metric, as described in this paper, and 
the point vortices for Diffvoi(IR.^) with an metric, and identify a number of questions that 
remain open. 



2 particle systems 

The Hamiltonian for the Euler equation is the kinetic energy 

-ff = - y {{u,m)) dxdy — 2 J ^^-^^^^^xdy. (5) 



Under the particle ansatz 

N 



m{x,t) ^^Pi{t)S{x ~ q,{t)), (6) 
{pi(t), qi{t) G for each i) the Hamiltonian becomes 

H = -Y.P.-pM\h-<i,\\)- (7) 



Solutions to Q of the form © obey Hamilton's equations for |7J, in which the components 
of Qi and Pi are canonically conjugate variables (see for further details). Here qi,. .. ,qN 
represent the positions of the TV particles, and pi, . . . ,pjv represent their momenta. 
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There are two immediate differences from the analogous system for point vortices in the 
Euler fluid equations: the number of degrees of freedom is 2N instead of iV, and, to get well- 
defined ODEs, the inertia operator A must be chosen so that G{Q) is finite. (For point vortices, 
G{r) = — log(r)/(27r), but the infinite self-energy of each particle can be ignored.) 

The equations of motion are 

TV 

i=i 

TV 

They have four conserved quantities: the energy H, the linear momentum X^i^i-Pi' ^'^d the 
angular momentum 1i ^ Pi- These are sufficient to ensure the integrability of the two 

particle problem, which is discussed in section ^ However, the dynamics of these two-body 
problems are very different from more familiar ones on flat configuration spaces, like the Kepler 
problem. Also note that Eqs. lO are not Galilean- invariant; their dynamics depends on the 
linear momentum. 



3 Blobs and kernels 

The Euler equation for Diff(M") and inertia operator A is given by Q and A related 
equation, the regularized Euler equation, is given by ^ together with u = * m (in place of 
where is the regularized Green's function 

G-,l,(x,x') — / G(x,x')'i/j{x') dx' . 



Typically, ipix) = is a 'blob', or approximate delta-function. 

For a given inertia operator A, the kernel G is the velocity corresponding to a blob /(|la;|| 
of momentum. It satisfies 

AG^fi\\x\\), 

which (assuming A is rotationally invariant) can be solved using the Hankel transform as 

u{x) — [ — — }-Jo{m\\x\\)m dm, 
Jo a(m) 

where 



/(™) = / f{r)Jo{mr)rdr 
Jo 

is the Hankel transform of /, and a(m) is the Fourier symbol of the operator A, e.g., d{m) = 
1 + for A^ 1- V^. 

Notice that considering a regularised Euler inertia operator with symbol d{m) together with 
blobs of symbol f{m) is identical to considering the standard Euler inertia operator with symbol 
d{m)/f(jn) together with delta-function particles. 

Consider first the metric associated with the inertia operator A2 — (1 — V^)^. Its kernel 
is G(r) = —rKi{r) /{An), where Ki{r) is a modified Bessel function of the first kind. Note that 
G(r) is negative for all r, reaching its global minimum (— l/(47r)) at r = 0. Thus, the velocity 
field induced by a single such particle is in the direction —p. We have found it more convenient 
to change the sign of all kernels, so that the induced velocity field is in the direction +p. Since 
the dynamics are reversible, this does not materially affect them. In this case, we work with 
G{r) = rKi{r)/{ATr). 
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Figure 1: Kernels Gfe(r) found in Section O for a family of iJ*^ metrics {k = 1,2,3,4, cx)) with 
delta-function momentum. 



For numerical simulations it is convenient to have an explicit form for G{r). In addition, 
it turns out that the dynamics of the 2-body problem depends sensitively on G{r) near r = 0. 
Therefore, we have worked with the family of H'' metrics with inertia operators Ak ■= (1 — 
Q,2y2-jfe g^jji-^ delta- function particles. Then the kernels Gk can be found explicitly in terms of 
modified Bessel functions. The first few are (in terms of the scaled length r :— r/a) 



27ra2Gi(r) 
2T:a^G2ir) 

and 27ra^G3(r) 



log(2/f) -7 + 0(7-2 \ogr), 



-rKi{r) 



2 

1 1 



8 



(21og(2/f) -27 + l)r 



;;2 



128 



(41og(2/r) + 5 - 47) f'' + O (r^) 



-firKoir) + 2Ki{f)) 
o 

1 1 1 

4 



f + ^ (4 log(2/f) -I- 3 - 47) f4 + O (r^) 

Choosing the length scale a ~ e/Vk we have formally, 

Ak Aoo := exp(-£2v2), 
with Fourier symbol a(m) = exp(— e^m^) and Green's function 



2^G^(r) = — exp(^- — 

That is, the Gaussian kernel is Green's function for an H°° metric, the limit of a family of H'' 
metrics. The family of kernels is shown in Figure Q 
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4 2 particle systems 



In this section we derive the reduced Hamihonian for the case where there are only 2 parti- 
cles. This enables us to describe the phase portraits, and to consider how they depend on the 
behaviour of the kernel, and hence on the metric. 



4.1 The reduced Hamiltonian 

Let = 2, so that the Hamiltonian (see Eq. (Tj) is 

H - pi-p2G{\\q^ - q2\\) + i(||pi|P + \\P2\?)G{0). 

We first change to the canonical Jacobi-Haretu centre of mass coordinates ^I] in order to reduce 
by the translational symmetry. Let 

d = pi+p2, P={pi-p2)/2, c={qi+q2)/2, Q = q^ - q2, 

or ^ ^ ^ ^ 

Pi = l^d+P, p2 = -d-P, q^=c+-Q, q2 = c--Q. 

Then 

H = - \\Pf) G{Q) + Qlldf + \\Pf^ G(0). 

Changing to symplectic polar coordinates with 

Q = {r cos6,r sm9), P = [pcosO — pg sin9/r,psm9 + pg cosO/r), 
where {r,p) and {0,pg) are conjugate variables, the reduced Hamiltonian becomes 

H = \\\df{G{Q) + G{r)) + [p' + ^) (G(0) - G(r)). (8) 

The conserved quantities are the momentum of the centre of mass, d, and the angular momentum 
pg —: /i. The phase portrait of this reduced Hamiltonian depends very sensitively on the 
behaviour of G near r = 0. 

First, consider the case when G is smooth at r = 0, so that G{r) — G(0) + 5G"(0)r'^ + 
^G^'^\Q)r^ + 0(r^). In this case, near r = 0, the Hamiltonian is 

(G(0)||d|p - G"[Q)^,^) + ^ (3G"(0)(||d|p - V) _ G(^)(0)^2^ + 0{r*). (9) 
Therefore, a separatrix intersects the p-axis at 

^ 4" " 12G"(0)'^ ■ 

(Typically, and in all the examples of this paper, G"(0) < and G^^\0) > 0.) The energy on 
the p-axis, and hence on the separatrix, is ^ (G(0)||(i|p — G"(0)^^). As r cx3, G(r) — > so 
H ~ G{0){p^ + -J II dp) and the separatrix has a horizontal asymptote at 

This can be clearly seen in Figure |21 where the phase portrait is given for the 2 body problem 
under the Aoo Gaussian kernel. The phase portraits of the 2 body problem under other metrics 
are given in Figures El and [3 and discussed in section ^21 

Thus there are three principal behaviours of pairs of particles: 



6 



123 + 5678 123 + 5678 
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d =0,(1 = 1 d =0,M =0 




Figure 2: Phase portrait for the 2 body problem with Gaussian kernel G{r) = exp(— r^/4). 
The Hamiltonian is given by Eq. Top: Nonzero linear momentum d. Bottom: zero linear 
momentum. Left: Nonzero angular momentum fi. Right: Zero angular momentum, for which 
the motion is confined to a line. 
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Scattering Both particles are free as t ^ ±00 
Capture Distance between particles r — > as i — > 00 
Ejection Distance between particles r ^ as t — > —00 

as well as degenerate capture/ejection orbits given by the separatrices. Inside the separatrix 
the two particles scatter, while outside the separatrix they either capture {pG"{0) < 0) or eject 
(pG"(0) > 0) each other. 

We now consider the criterion for a capture to take place. Equation H10() tells us that two 

particles initially very far apart will capture each other if > -jUdlP — ^g(oj M^- This criterion, 
depending on the three parameters p, \\d\\, and /i, can be expressed more intuitively as follows. 
We introduce the angle of approach ip := Zpip2, the speed ratio Ibi||/|b2||, and the separation 
distance D. D is the distance of closest approach of two particles initially at qi^2 moving in 
straight lines with direction pi^2- Clearly, as D increases then the particles cease to interact, so 
we will seek the critical value of D for a capture to occur. 

We have D = \\Q\\ sin(ZQP) and fi = Q x P = \\Q\\\\P\\ sin(ZgP) = D\\P\\. As r ^ cx), 
p^ ||-P|P — jWpi — so the capture criterion becomes 

-PI ■P2 = -Ibl|l|b2|| COSV- > -^^D^P\\\ 

or, using the law of cosines ||2P|p = — = + IIP2IP — 2||pi|| ||p2|| cost/', 

cos^<^fM + « _1 5<0, 
l + 2i^V 11^211 bill/ 4 



where 



8G(0) 



That is, a separation distance D of less than —2G{0)/G"{0) is necessary for a capture to be 
possible. Particles with D = (equivalently, 11 = 0) are captured when IV'I > 7r/2 (a 'head-on' 
approach) and scattered when < 7r/2 (a 'glancing' approach), for all Ibi||/|b2||; increasing 
D and ||pi|/||p2|| moving away from 1 both restrict captures to more nearly head-on approaches. 
Figure 121 shows the minimum required approach angle "ip for each D. 

During a capture the particles get infinitely close, r ^ 0. Although p —00, the actual 
speed of approach of the particles, r, tends to 0. The particles do not collide in a finite time but 
rather become bound together, or stick to each other. For large p and small r, we can use the 
approximate form of the Hamiltonian for r « to get the approximate equations of motion: 

r = -G" {0)pr'^ , p^G"{0)p^r, 

with solution 

r(i)=roe-"*, p{t) ^ poe"\ (11) 

where the constant a — G" (O)poro, corresponding to a capturing orbit for pqG"{0) < and an 
ejecting orbit for poG"{0) > 0. We call a pair of particles in the limiting state of Hll|l a dipole, 
and explore them further in section [S| 
For a scattering orbit, using 



dH _ G{0)-G{r) 
ope r'^ 



we see that the scattering angle 9\°^^ ranges over (0,oo) as the minimum separation of the 
particles decreases from 00 to 0. (As r 0, = 0(1) and p — 0(r), so the phase 9 accumulates 
without bound.) Again, because close particles move very slowly, there can be an arbitrarily 
long delay during the close approach — the particles appear to stick together for a while. 
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0.2 0.4 0.6 O.a 1 1.2 1.4 1.6 1 .S 2 

D 

Figure 3: Minimum required approach angle ip G (f-jTr) for a capture to take place for a pair 
of initially distant particles with speeds ||pi,2|| and separation distance D. Here we have taken 
G(r) — e^^ for which the maximum possible separation distance in a capture is 2. That 
this equals the range of G is something of a coincidence, because the critical D depends only 
on G(0) and G"(0). 



4.2 Particle dynamics under different metrics 

We can consider different kernels (i.e., different functions G{r)) and use them to discuss the 
particle dynamics for each metric in turn. In all cases, we consider the Helmholtz-style inertia 
operator described previously: Ak (1 — a^V^)*''. We begin with the metric, and progress 
through the different H'^ metrics. 

Delta-function particles for the inertia operator Ai {H^ metric) are not well-defined because 
they do not have a finite self-induced velocity (the same situation occurs in the 3D Euler fluid 
equations when the vorticity is concentrated on a curved filament). To use the delta- function 
particles, it is necessary to regularize the equations. One approach to this would be to simply 
set G(0) — 0, so that the self-induced velocity of each particle is 0. We consider this below 
(section I4.3|l , but we first consider the standard regularization of introducing a smooth blob 
function /(x), so that we consider the operator A = Ai/ f{x), as we described earlier. 

As mentioned in the introduction, we consider the Gaussian blob f{x) = -^e"^ /e^, which 

tends to 5(x) as £ ^ 0. In this case, j{m) — -^e^^ ™ The kernel for this blob is smooth 
with the inertia operator = 1 — , so we apply the results of section 14.11 and expand for 
small X as: 

2-G(0 = ^ + + i^) + - ^ + 1,C) x^ + 0ix% 

where C — ^ exp(i£^) Eii(ie^). 

As e we have G(0) - (-21oge - 7 + log4)/(47r), G"(0) l/(27re2), and G""(0) - 

3/(27re'*), and so the separatrix intersects r — at p = 0{^/e), with a horizontal asymptote 
at p = 0(/i/(e(log |e|)^/^). As e — > 0, providing that /i 7^ 0, the separatrix moves outward and 
particle captures become rarer. 

The metric is smooth enough for delta-function particles to be well-defined, with a finite 
self-induced velocity. The phase portrait for this case is shown in Figure. 01 However, G2(f) is 
not twice differentiable at r = and so the above analysis of the phase portrait of the two body 
problem is not valid. In fact, the singularity of G2{r) is sufficient to prevent the occurrence of 
capture orbits; when ^ all orbits are scattering (see Figure 0J. However, two particles with 
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d = 1,M = 1 d = 1,M=0 




123 + 5S7S 1234567a 



d =0, M = 1 d =0, M =0 




r r 



Figure 4: Phase portrait for the 2 body problem with inertia operator A2 = [l — V^)^. The 
Hamiltonian is given in Eq. JHJ. Unless the angular momentum /i is zero (in which case the 
motion is confined to a line), all orbits are scattering. 



large initial relative momentum p can approach arbitrarily closely, and achieve arbitrarily high 
p during the close encounter (which makes them behave a lot like capturing orbits); but they 
do eventually separate. A typical such encounter is shown in Figure 

The metric is smoother, but G3(r) is still not 4 times differentiable at r = 0, so the above 
separatrix analysis is still not valid. However, one can check that the 2 body problem does now 
have a separatrix that divides capture and scattering orbits, although it is not a smooth curve 
at r = 0. 

For H'' [k > 4) and H°°, the separatrix is as described above. Figure [3 shows the phase 
portrait for the metric. A typical scattering orbit is shown in Figure O in contrast to the 
case, the relative momentum 2P — pi ^ P2 is bounded over all scattering orbits. 
In no case are there any periodic orbits in the two body problem. 
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Figure 6: A typical scattering orbit for the 2 body problem with inertia operator Aoo- The 
relative momentum 2P = p\—p2 remains bounded during the close approach. 
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4.3 An alternative regularization for 

As was mentioned above, there is another possible regularization of the A^-particle problem 
with the metric, which is simply to set the self-induced velocity of each particle to zero, 
i.e., G'(O) = 0. This seems drastic, when the 'correct' self- induced velocity of a delta- function is 
infinite, but it still corresponds to a consistent discretization of the PDE I^J in the limit of a 
large number of particles spaced over a curve or an area. (The contribution to the velocity at 
a point from nearby momentum p is p{r)Gi{r)r dr for a curve, and J^J^^ p{r,9)Gi{r)r^ dr 
for an area; in both cases the logarithmic singularity in G'i(r) is weak enough that the nearby 
momentum does not contribute.) 

Setting G(0) = changes the character of the two body problem completely. It now features 
periodic orbits and scattering orbits, but no capturing orbits (see FigureEj). However, numerical 
simulations of the three body problem in this case indicates that orbits can reach a singularity 
in a finite time, via a mechanism in which 2 particles orbit closer and closer, under regulation 
from the approach of the third particle. 



5 Capture orbits and dipoles 

Consider the A^-body problem together with a smooth kernel so that 2-body captures are pos- 
sible. Suppose that particles 1 and 2 become close enough to capture each other and that there 
are no other nearby particles. We will call the limiting state of such a captured pair of particles 
a 'dipole'. In centre-of-mass coordinates for particles 1 and 2 only, the dipole is described by its 
centre of mass c, its momentum d, and its orientation Q (equivalently by (r, 9)) with conjugate 
momentum P (equivalently (ji,pg)). From Eq. Elwe get that the Hamiltonian of a lone dipole 
is given by 

Hi,2 :=-iG"(0)(rV+P^) + ^G(0)||d|r. 

The motion of a lone dipole is given by 

c{t)^CQ + tG{0)do, d{t)=do, r{t)^roe-°'\ p{t) ^ poe'^K, 0{t) = 9o + ojt, (12) 

where the frequency lu = -G"{0)fi and a = G"{0)poro > (G"(0) < 0, po < 0). However, 
other particles couple to all the degrees of freedom of the dipole so that it continues to undergo 
internal evolution. 

Writing Hi jy for the Hamiltonian of particles 1 through N, the terms in Hi jy coupling 
the dipole to particle j are: 



Hf := p,- 



Id + P\G (\\c - q, + \q\\\ + (\d -p\g(\\c- q, - \q\\ 



[d ■ p,)G{\\c ~ II) -I- {p, ■ P){(c - q,) ■ 0)-^"^ 



plus terms of order 0(6""*), which we omit if the dipole is taken to be in its limiting state. The 
first term in H'^ is the Hamiltonian for a single particle at position c with momentum d, while 

the second (dipole) term generates an 0(1) contribution to all of d, pj, qj, Q/Q, and P/P. 
Identifying qN+i = c and pn+i = d, the Hamiltonian for a dipole and N — 2 singleton particles 
is 



11 1 ^ 

2^1,2 + 2^3,...,Ar+i + 2 H 



i=3 
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Figure 7: Phase portrait for the 2 body problem with inertia operator Ai and regularization 
G(0) = 0. The dynamics is somewhat reminiscent of traditional 2-body problems like the Kepler 
problem of Newtonian gravity: there are no capture orbits; for nonzero angular momentum 
(left) there are scattering and periodic orbits, while for zero angular momentum (right) genuine 
collisions are possible. 
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This suggests that the overall evolution of N particles with a smooth kernel is described by 
(i) free, straight-line motion of all particles when all particles are well-separated; (ii) scattering 
and interactions when distances are moderate; and (iii) capture or sticking together of 2 particles 
when they come within range, after which they continue to evolve as a dipole (and interact with 
the other particles). It is also possible for more than 2 particles to simultaneously capture each 
other. We have not analyzed this situation but it seems to follow the same course as the 2 body 
capture, with all inter-particle distances tending to exponentially fast. It may also be possible 
for the influence of the other particles to cause a dipole to separate, although we never observed 
this happening. 

If the boundary conditions are taken to be periodic, this indicates that most orbits of the 
iV-particle problem end up with all particles stuck together in one lump. 

A typical orbit for 20 particles in the plane is shown in Figure |S1 The initial positions 
and velocities are (normally) randomly distributed. Here the particles have been captured into 
groups of 9, 3 (with a 4th about to join), 2 (with 2 more possible captures), and 3 singletons. 

We now analyze the influence of a dipole on a third particle with zero momentum — a 'test 
particle'. The test particle does not affect the motion of the dipole. We only have to flnd the 
position q of the test particle, which obeys 

q = d{t)G{p) + P{{c{t) - q) ■ Q)fip), (13) 

where c{t) and d{t) are given by Eq. (|12|l . p = \\c — q\\ is the distance from the dipole to the 
test particle, and f{p) = G'{p)/p. 

One solution of H13() is q{t) — c{t), in which the test particle sits on top of the dipole. We 
shall see that dipoles can capture test particles into this state. 

First consider the case c? = of a motionless dipole located at c = 0. We compute 

Poro cos ut — p sinuit 
Poro sin tot + p cos uit 

- sin/? 
cos/? 

r cos ujt 
r sin Lut 

where = p^r^ + p^ and (3 = cut + (p, and the phase (p is defined by cos0 = p/A, s'mc 
—Pofo/A. Eq. l(T^ becomes 



rP = 
^A 



q = —A (qi cos Lut + q2 sin Lot) f{p) 



— sin (3 

cos (3 



note the e"* factors cancel out and the interaction is characterized by the parameters ropo a-nd 
p of the dipole. The form of these equations suggests going into a rotating frame by 

cos f3 sin /3 



y I sin/3 cos/?/"^' 
which leads to a drastic simplification to the planar, autonomous system 

yi y2, ^^^^ 

h = -{mn + roPoy2)f{p) - ujyi- 

Note p = \\q\\ — \\y\\. The phase portrait of this system is shown in Figure [Till 

We now investigate the stability of the fixed point y = under the assumption that G is 
smooth, so that f{p) = G"(0) + ^G^^\0)p^ + ... and ((HI) becomes 

m = ojy2 = -G"{0)py2 

2/2 = roPoG"{0)y2 - Ig^'H^M + yl) + oiWyf). 
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Figure 8: Typical orbit for 20 randomly placed particles under a Gaussian kernel. (Top: par- 
ticle paths, with final particle positions indicated by a dot; bottom: snapshots of the particle 
positions at 16 different times.) 
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Figure 9: Attraction of a test particle to a dipole at the origin. (The particle moves (under Eq. 
(|12|) in the direction of the rotating direction vector of the dipole, with the cusps corresponding 
to the dipole pointing directly at the particle.) 
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The eigenvalues of the fixed point y = are (eigenvector (1,0)"^) and — po'"oG"(0) < 
(eigenvector (/i,po^o)"'")- The fixed point is nonhyperbolic, but a routine application of cen- 
tre manifold theory then gives a centre manifold located at j/2 — JVi + 0{yf) with 7 = 
— /iG'^^-'(0)/(6po^oG"'(0)), and, on the centre manifold, the dynamics is governed by the reduced 
equation 

JC 3 r 

yi^oy-^, d =u;j = - 

oporo 

Recalling that po < for a dipole, we have that ?/ = is asymptotically stable if G("''(0) > 
and unstable if G'-'*-' (0) < 0. The fixed point is stable for the Gaussian kernel G(p) — e^^ . 

That is, the dipole for a Gaussian kernel attracts nearby test particles, albeit very slowly, 
since yi{t) = (yi(0)~^ — 25t)~^^'^. A typical test particle far from the dipole moves initially 
slowly (since the velocity field is exponentially small there), then quickly falls onto the centre 
manifold, then orbits the dipole at the same frequency as the dipole itself, maintaining a phase 
shift of (j), at a gradually diminishing distance. See Figures El and EH 

Finally, we consider a moving dipole, which we take (e.g., by rotating space and by scaling 
time) to be moving along the x-axis at unit speed, c{t) = t{l, 0). A test particle with sufficiently 
large |92(0)| (i.e., with the test particle far enough from the path of the dipole) will not be 
entrained as the dipole passes. However, the critical value of (72(0) will also depend on qi{0), 
which determines the phase of the dipole as it approaches. A numerical simulation (see Figure 

reveals a surprising asymmetry where, for an anticlockwise (^ > 0) dipole, the entrainment 
is independent of the phase for 92(0) < and strongly dependent on the phase for 92(0) > 0. 
All particles out to (72 = —1.536, on the right of the dipole, are entrained, but on the left the 
critical initial distance varies from 0.6 to 1.15. 



6 Discussion 

The particle dynamics in Diff(M^) with a H'' metric {k > 2) that we have considered in this 
paper show striking contrasts to the more familiar dynamics of point vortices for Diffvoi(lR^) 
with an metric. Of particular note are: 

• our particles are vector particles, point vortices are point particles; 

• for our particles, the 2 body problem is integrable, for point vortices the 3 body problem 
is integrable; 

• our particles exhibit short-range interactions u{r) ~ g~'" 7.^^1/2 (^j. _j. qq^^ point vortices 
exhibit long-range interactions u{r) ^ 1/r; 

• our particles have no periodic 2-body orbits, whereas most 2-body orbits are periodic for 
point vortices; and 
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Figure 10: Phase portrait of Eq. H14(l for a test particle in the field of a motionless dipole in 
rotating coordinates. Here, /i = 1, ropo = ~li Sind G{p) — e~P . 

• the 2-body dynamics of Diff(R^) depend sensitively on the metric, whereas for point vor- 
tices, the qualitative 2-body dynamics are independent of the metric. 

Clearly, we have only scratched the surface of the rich dynamics of the system (0 . Further 
dynamical questions to be considered include 

• a full study of the 3 body problem; 

• determination of the long-time existence of solutions; 

• a classification of the limiting states of the A^-body problem; 

• a study of the limiting (W-pole') state of the capture of N particles; and 

• a consideration of the geometry of the Riemannian manifold whose geodesies are governed 
by ©, 

while broader questions include 

• the significance of the sensitivity to the metric for image matching applications; and 

• the convergence of the point particle approximation to solutions of the PDE H1I2() . 
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Figure 11: Entrainment of a field of test particles by a moving dipole. The dipole (whose 
direction is indicated by a stick) is moving to the right with speed 1. Particles shown as small 
dots are entrained by the dipole, while particles shown as larger dots are left behind. 
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